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Abstract 

Using full-potential, density-functional-based methods we have studied electron-phonon interac- 
tion in ZrB2 and TaBi in P6/mmm crystal structure. Our results for phonon density of states 
and Eliashberg function show that the electron-phonon coupling in ZrB2 is much weaker than in 
TaBi- In particular, we find that the average electron-phonon coupling constant A is equal to 0.15 
for ZrB2 and 0.73 for TaBi- The solutions of the isotropic Eliashberg gap equation indicate no 
superconductivity for ZrB2 but a superconducting transition temperature T c of around 12 K for 
TaB2 with //* ~ 0.16. 
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I. INTRODUCTION 



The search for superconductivity in the transition-metal diborides ZrB^ and TaB 2 has 
proven to be elusive since the first reports of superconductivity in ZrB 2 



Bo and TaB 2 has 
Q and TaB 2 Q 



with superconducting transition temperatures T c = 5.5 K and T c = 9.5 K, respectively. The 
subsequent experimental efforts jj], U] to make ZrB 2 and TaB 2 superconducting by applying 
external pressure and/or hole doping have also proven to be futile. Such an outcome for 

ZrB 2 and TaB 2 seems somewhat surprising given that under similar conditions Nb x B 2 

an 
and Moo.96Zr OM B 2 |5j] are found to superconduct with T c = 8.5 K and Tc = 5.9 K, 

respectively. 

To be able to understand the superconducting properties or lack thereof of the transition- 
metal diborides ZrB 2 and TaB 2 within the conventional BCS theory of superconductivity 
one needs to know the electron-phonon interaction in these systems. In particular, the 
strength of the electron-phonon interaction in ZrB 2 and TaB 2 , as reflected in the average 
electron-phonon coupling constant A, can be used to predict the possibility of superconduc- 
tivity in the diborides. Experiments, based on point-contact spectroscopy [f3,Ql, indicate A 
to be less than 0.1 for both ZrB 2 and TaB 2 , and thus preclude the possibility of supercon- 
ductivity in these diborides. Here, we like to point out that the experimental determination 
of A, as carried out in Refs. jfj Q], is more reliable for ZrB 2 than for TaB 2 . 

Previous theoretical estimates of A Q, 0], based on only zone-center sampling of the 
electron-phonon matrix elements, also indicate ZrB 2 and TaB 2 to have a weak electron- 
phonon coupling with A < 0.2. However, it turns out that in the transition-metal diborides 
the estimates of A based on zone-center sampling of the electron-phonon matrix elements 
corresponding to the optical E 2g mode is not reliable P, Thus, a detailed analysis 

of the electron-phonon interaction using a more representative Brillouin zone sampling of 
the electron-phonon matrix elements is clearly needed to better understand the lack of 
superconductivity in ZrB 2 and TaB 2 . The present work is a step in that direction. 

Using density-functional-based methods we have studied (i) the electronic structure, (ii) 
the phonon density of states (DOS), (iii) the electron-phonon interaction, and (iv) the so- 
lutions of the isotropic Eliashberg gap equation for ZrB 2 and TaB 2 in PQ/mmm crystal 
structure. We have calculated the electronic structure of ZrB 2 and TaB 2 with experimental 
lattice constants a and c, as given in Table I, using full-potential linear muffin-tin 
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orbital (LMTO) method. For studying the electron-phonon interaction we used the full- 
potential linear-response program of Savrasov [ll|, [l^] to calculate the dynamical matrices 
and the Hopfield parameters, which were then used to calculate the phonon DOS, F(lj), the 
electron-phonon coupling A, including the partial A q , and the Eliashberg function, a 2 F(u), 
for ZrE>2 and TaB?. Subsequently, we have numerically solved the isotropic Eliashberg gap 
equation Q,Q,Q for a range of /x* to obtain the corresponding superconducting transition 
temperature T c . 

Based on our calculations, described below, we find that the electron-phonon coupling in 
ZrB 2 is much weaker than in TaB 2 . In particular, we find that the average electron-phonon 
coupling constant A is equal to 0.15 for ZrB 2 and 0.73 for TaB 2 . The solutions of the 
isotropic Eliashberg gap equation indicate no superconductivity for ZrB 2 but a supercon- 
ducting transition temperature T c of around 12 K for TaB 2 with fi* ~ 0.16. 



II. THEORETICAL APPROACH AND COMPUTATIONAL DETAILS 

Within Migdal-Eliashberg theory of superconductivity , first-principles calcula- 

tions of superconducting properties require the knowledge of (i) the ground-state electronic 
structure, (ii) the vibrational spectrum and (iii) the electron-phonon matrix elements of 
the solid. The phonon spectrum and the electron-phonon matrix elements are then used to 
calculate the Eliashberg functions from which the superconducting properties of the mate- 
rials can be calculated. In particular, by solving the fully anisotropic or the isotropic gap 
equation the superconducting transition temperature can be calculated. 

The density-functional theory provides a reliable framework for implementing from first- 
principles the Migdal-Eliashberg approach for calculating the superconducting properties of 
metals as outlined above. Such an approach has been implemented by Savrasov using the 
jMTO formalism and the linear response method, the details of which are given in Refs. 

. In the following, we briefly outline the approach used in the present calculations. 
We follow the notation used in Refs. [n|, Q|- 

We calculate the ground-state electronic structure of the transition-metal diborides ZrB 2 
and TaB 2 using the full-potential linear muffin-tin orbital method in the generalized-gradient 
approximation of the density-functional theory. The linear response approach is used to eval- 
uate the dynamical matrices and the electron-phonon matrix elements within the density- 
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functional theory as implemented using the LMTO formalism. It amounts to self-consistently 
evaluating the changes in the electronic structure due to atomic displacements associated 
with the phonon mode a> q „. Once the changes in the wave function, density and the effective 
potential are known, the phonon density of states 

FM = (1) 
and the electron-phonon matrix elements 

C^=<k + qj-'|^K//|kj>, (2) 

can be evaluated. The electron-phonon matrix elements can be interpreted as the scattering 
of electron in state |kj > to state |k + qj' > due to perturbing potential 5 qu V e ff, which arises 
due to the phonon mode ofy . Note that Eq. [2] has to be corrected for the incomplete basis 
set as described in Ref. From the electron-phonon matrix elements we can calculate 

the phonon linewidth 7 qi , 



7 qv = 27Ru qi/ ^ |£k+q;',kj| 2( K e ki - e F )5(e k+qi * - e F ). (3) 

Occasionally, it is useful to define a phonon mode dependent electron-phonon coupling A qi/ 
as 



where N(sp) is the electronic density of states at the Fermi energy. Now we can combine 
the electronic density of states, phonon spectrum and the electron-phonon matrix elements 
to obtain the Eliashberg function a 2 F{uj) defined as 

27r7V(e F ) q „ 

Finally, the Eliashberg functions can be used to solve the isotropic gap equation [l3j, LLJ, Il5| 



\uj n ,\<w c 

A(iu n )= fn'S(n,n')A(iuj n/ ) (6) 

n' 

to obtain the superconducting properties such as the superconducting transition temperature 
T c . The function S(n, n') used in the gap equation is defined by 

S(n, n) = \{n - n) - /i* - 8 nn > s n s n "K n ~ n ") ( 7 ) 



Table I: The experimental lattice constants for ZrB 2 and TaB 2 used in the calculations. 



nil cwj 

allU V 


n la 11 1 

\Jj Id, LI . J 


on 


ZrB 2 


5.992 


1.1142 


TaB 2 


5.826 


1.0522 



and f n = l/\2n + 1| with s n representing the sign of u n . The electron-phonon coupling A(^) 
is given by 

roo Or,} 

X(y) = / dua 2 F(u) — - (8) 

JO UJf, + UJ Z 

Before describing our results in detail, we provide some of the computational details 
which are similar to our study of the electron-phonon interaction in MgB 2 , NbB 2 and TaB 2 

The charge self-consistent full-potential linear muffin-tin orbital calculations were carried 
out with the generalized gradient approximation for exchange-correlation of Perdew et al 
Hi and 484 k-points in the irreducible wedge of the Brillouin zone. For TaB 2l the basis 
set used consisted of s, p, d and / orbitals at the Ta site and s, p and d orbitals at the B 
site. In the case of ZrB 2 , we included s, p and d orbitals at the Zr site. In all cases the 
potential and the wave function were expanded up to l max = 6. The muffin-tin radii for Ta, 
B, and Zr were taken to be 2.5, 1.66, and 2.3 atomic units, respectively. 

The calculation of dynamical matrices and the Hopfield parameters for ZrB 2 and TaB 2 
were carried out using a 6 x 6 x 6 grid resulting in 28 irreducible q-points. The Brillouin 
zone integrations for charge self-consistency during linear response calculations were carried 
out using al2xl2xl2 grid of k-points. The Fermi surface was sampled more accurately 
with a 36 x 36 x 36 grid of k-points using the double grid technique as outlined in Ref. 
[l^ . In case of TaB 2 , for a couple of q-points we had to sample the Fermi surface using a 
48 x 48 x 48 grid of k-points to get converged acoustic-mode frequencies. 



III. RESULTS AND DISCUSSION 



In this section we describe the results of our calculations of the electronic structure, the 
linear response and the solutions of the isotropic Eliashberg gap equation for ZrB 2 and 
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TaB 2 . Our results for ZrB 2 and TaB 2 are described in terms of (i) the density of states 
(DOS), (ii) the phonon density of states F(uj), (iii) the Eliashberg function a 2 F(u), and (iv) 
the superconducting transition temperature T c obtained from the solutions of the isotropic 
Eliashberg gap equation. 



A. The Density of States 

The electronic structure of ZrB 2 [3] and TaB 2 jj, Q, Q| has been studied previously. 
Our full-potential results for the electronic structure of ZrB 2 and TaB 2 are in agreement 
with earlier calculations. For example, in Fig. 1 we show the total densities of states for 
ZrB 2 and TaB 2 calculated at the experimental lattice constants using the full-potential 
LMTO method as described earlier. The agreement with earlier calculations of the densities 
of States for ZrB, } and TaB 2 Q „ exceHen, We fa d tt at the fota! DOS at the Fe rmi 
energy is equal to 3.72 st/Ry and 12.52 st/Ry for Zr_B 2 and TaB 2 , respectively. In ZrB 2 , 
the relatively low DOS at the Fermi energy leads to a weak electron-phonon coupling in this 
system. However, at the Fermi energy the d-electrons are present in substantial amount in 
both ZrB 2 and TaB 2 , indicating a more active role for Zr and Ta atoms in determining 
the lattice- dynamical as well as the possible superconducting properties of these diborides. 



B. The Phonon Density of States 

In Fig. 2 we show the phonon DOS F(u) of ZrB 2 and TaB 2 calculated using the full- 
potential linear response program as described earlier. As is clear from Fig. 2, the phonon 
DOS for both ZrB 2 and TaB 2 can be separated into three distinct regions. Based on our 
analysis of the eigenvectors, we find that the first region, with a peak in phonon DOS at 
32meV in ZrB 2 and 22meV^ in TaB 2 , is dominated by the motion of the transition-metal 
atoms Zr and Ta respectively. Note that the shift of the first region in the phonon DOS 
towards lower frequency for TaB 2 in comparison to ZrB 2 is due to the higher mass of Ta. 
In the second region, the phonon DOS around 60 — 70 meV in ZrB 2 (TaB 2 ) results from the 
coupled motion of Zr (Ta) and the two B atoms. However, the peaks in the phonon DOS 
at 78 meV in ZrB 2 and 75 meV in TaB 2 correspond to the in-plane B — B motion. The 
phonon DOS in the third region, which extends from 94meV to 102 meV in ZrB 2 and from 
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Figure 1: The total density of states of ZrB 2 and TaB 2 calculated using the full-potential linear 
muffin-tin orbital method as described in the text. The vertical dashed line indicates the Fermi 
energy. 
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Figure 2: The phonon density of states F(u) of ZrB 2 and TaB2 calculated using the full-potential 
linear response method as described in the text. 



100 me V to 107 meV in TaB 2 , results from the displacements of all the three atoms. Not 
surprisingly, the phonon DOS of ZrB 2 and TaB 2 are similar to the phonon DOS of NbB 2 
if we make allowance for the difference in the masses of the transition-metal atoms |lfl| . 



7 



0.4 
0.3 
0.2 
? 0.1 

P 

8 1.2 
0.8 
0.4 
0, 



t 1 1 r 




v I I 




20 40 60 
co(meV) 



100 120 



Figure 3: The Eliashberg function a 2 F(u>) of ZrB 2 and TaB 2 calculated using the full-potential 
linear response method as described in the text. 



Table II: The calculated Hopfield parameter 77, the average electron-phonon coupling constant A, 
the root mean square u> rms and the logarithmically averaged u>i n phonon frequencies for ZrB2 and 
TaB 2 . 



alloy 


77 (mRy / a.u. 2 ) 


A 


Urms (K) 


W/ n (K) 


ZrB 2 


76 


0.15 


734 


585 


TaB 2 


279 


0.73 


593 


356 



C. The Eliashberg Function 

The main purpose of the present study is to examine the strengths of the electron-phonon 
interaction in ZrB 2 and TaB 2 in order to better understand the lack of superconductivity 
in these diborides. To this end, we show in Fig. 3 the Eliashberg function a 2 F{uj) of ZrB 2 
and TaB 2 calculated as described earlier. A comparison of a 2 F{ui) of ZrB 2 and TaB 2 shows 
that the electron-phonon coupling in ZrB 2 is much weaker than in TaB 2 . In particular, the 
average electron-phonon coupling constant A is equal to 0.15 for ZrB 2 and 0.73 for TaB 2 . 

From Fig. 3, we also see that in the case of ZrB 2 the phonon mode with peak at 78 meV, 
which corresponds to the in-plane B — B motion, couples to the electrons more than the 
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Figure 4: The superconducting transition temperature T c as a function of fi* for TaB 2 as obtained 
from the isotropic Eliashberg gap equation. 

phonon mode with peak at 32 meV. However, in TaB 2 it is the phonon mode corresponding to 
the displacements of the transition-metal atom Ta, with a peak at 22 meV, which couples 
more strongly to the electrons. Such a change in a 2 F(u) clearly indicates a more active 
role for the transition-metal atoms in deciding the normal as well as the superconducting 
state (if any) properties of the transition-metal diborides. Our calculated a 2 F(u) for ZrB 2 
agrees with the point-contact spectroscopy results of Refs. [jj. However, in our opinion, 
the experimentally obtained a 2 F(u) for TaB 2 in Ref. (3| is underestimated, and it is in 
disagreement with the present result. In Table II we have listed the Hopfield parameter 77, 
the electron-phonon coupling constant A, and the various averages of the phonon frequencies 
for ZrB 2 and TaB 2 . 



D. The Superconducting Transition Temperature 



The possibility of superconductivity in ZrB 2 and TaB 2 within the present approach can 
be checked by solving numerically the isotropic gap equation Q| using the calculated 
Eliashberg function a 2 F{uj). The results of such a calculation for TaB 2 are shown in Fig. 4 
for a range of values of /x*. From Fig. 4 we find that for /i* ~ 0.16 the T c for TaB 2 is equal 
to ~ 12 K. A similar calculation for ZrB 2 using the calculated Eliashberg function yields a 



9 




u 20 40 60 
0)(meV) 



100 120 



Figure 5: The phonon density of states F(u>) of ZrB 2 and TaB 2 obtained with 12 q points (dotted 
line) and 28 q points (solid line) using the full-potential linear response method as described in the 
text. 

T c of less than 10~ 5 K. 

Our results show that ZrB 2 is unlikely to show superconductivity, in agreement with ex- 
periments I20J. In addition, the calculated Eliashberg function for ZrB 2 agrees well with the 

no 

corresponding Eliashberg function obtained from point-contact spectroscopy p,l2fl- However, 
our calculations predict TaB 2 to superconduct with a T c ~ 12 K but so far no superconduc- 
tivity Q has been found (except for Ref. Although our calculated values of a 2 F(u) 
and consequently A can be improved upon by including more q-points in the linear response 
calculations but it is unlikely to change the main conclusions of the present work. Thus, the 
reasons for the lack of superconductivity in TaB 2 need to be explored further. 



E. Convergence of Phonon Density of States, Eliashberg Function and Supercon- 
ducting Transition Temperature 

The first-principles study of superconducting properties of solids as carried out above is 
computationally demanding. In particular, the requirement of electronic self-consistency for 
phonon wave vector q and mode v makes the computational effort prohibitive unless one 
makes a judicious choice of q vectors. However, one must ensure that the calculated results 
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Figure 6: The Eliashberg function a 2 F(co) of ZrB 2 and TaB 2 obtained with 12 q points (dotted 
line) and 28 q points (solid line) using the full-potential linear response method as described in the 
text. 

are converged with respect to the number of q vectors. To demonstrate the convergence of 
our calculated phonon density of states F(u), Eliashberg functions a 2 F{uj) and the super- 
conducting transition temperature T c , we have carried out two sets of calculations using the 
double-grid technique as outlined in Ref. ll|, We have used two q grids (i) a 4 x 4 x 4 
grid with 12 irreducible q-points and (ii) a 6 x 6 x 6 grid with 28 irreducible q-points. All 
the results described so far correspond to the second grid with 28 q-points. Below we show a 
comparison of F(uj), a 2 F{oj) and T c calculated using the two grids with 12 and 28 q-points, 
respectively. 

In Fig. 5 we compare F(u) calculated using the two grids with 12 and 28 q-points, 
respectively. In general, dynamical matrices and consequently the phonon spectrum do not 
require a very fine k-mesh for convergence. However, the increase in the number of q-points 
has understandably reduced the peak heights at u approximately equal to 30 meV , 75 meV 
and 95meV for ZrB 2 and 20meV, 70meV and 100 meV for TaB 2 . 

The convergence of a 2 F{ui) is shown in Fig. 6, where we have plotted a 2 F{ui) calculated 
using the two grids with 12 and 28 q-points, respectively. For ZrB 2 the differences are small 
in magnitude except at around 20 meV, while for TaB 2 the effects due to small number of 
q-points are reduced with the use of the 28 q-point grid. We also see that the integrated 
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Figure 7: The superconducting transition temperature T c as a function of /x* for TaB 2 as obtained 
from the isotropic Eliashberg gap equation. The dashed and the solid lines correspond to the 
calculations done with 12 and 28 q points, respectively. 

values such as the electron-phonon coupling constant A do not change appreciably between 
the two sets of calculations (A changes by ~ 0.01 for ZrB 2 and ~ 0.005 for TaB 2 ). 

As F(u) and a 2 F{uj) do not change significantly between the two sets of calculations, 
we do not expect a significant change in the solutions of the isotropic gap equation. Con- 
sequently, the calculated T c should be similar in the two cases. In Fig. 7, we show the T c 
obtained for the two sets of calculations, and the difference of less than 1° K between the 
two sets of calculations confirms the convergence of F(cu), a 2 F(u) and T c . 

IV. CONCLUSIONS 

We have studied electron-phonon interaction in ZrB 2 and TaB 2 in PQ/mmm crystal 
structure using full-potential, density-functional-based methods, we find that the electron- 
phonon coupling in ZrB 2 is much weaker than in TaB 2 . In particular, we find that the 
average electron-phonon coupling constant A is equal to 0.15 for ZrB 2 and 0.73 for TaB 2 . 
The solutions of the isotropic Eliashberg gap equation indicate no superconductivity for 
ZrB 2 but a superconducting transition temperature T c of around 12 K for TaB 2 with fj,* ~ 
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0.16. 
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